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Abstract 

We present the exact solution of the Osp{l\2) invariant magnet by the Bethe 
ansatz approach. The associated Bethe ansatz equation exhibit a new feature by 
presenting an exphcit and distinct phase behaviour in even and odd sectors of the 
theory. The ground state, the low-lying excitations and the critical properties are 
discussed by exploiting the Bethe ansatz solution. 
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The concept of integrability based on the Yang-Baxter relation can be extended in 
order to include in its algebraic structure variables possessing both commuting and an- 
ticommuting rules of permutation. This makes it possible to construct integrable lattice 
models in which their bond configurations admit bosonic and fermionic degrees of free- 
dom This notion is then systematized by relating the Boltzmann weights of the lattice 
model to the i?-matrix solutions of a graded version |T]] of the Yang-Baxter equation. This 
generalization is necessary in order to take into account the Grassmann parities coming 
from the interchange of fermions of the theory. Moreover, for a typical system consisting 
of n species of bosons and m species of fermions, their solutions are believed to be found 
as invariants under the Sl{n\m) and Osp{n\2m) superalgebras 0. 

The exact solution of the Sl{n\m) invariant systems has been studied in many different 
contexts in the literature. For instance, the Sl{n\m) property is resembled much on the 
Perk-Schultz vertex model ^ which possesses an extra (discrete) parameter playing 
the role of the Sl{n\m) Grasmann parities. More recently, these solutions have also been 
examined in strongly correlated eletronic systems, being the one dimensional supersym- 
metric t — J model the typical example [0, 0]. However, for the Osp{n\2m) series, not 
much is known concerning its Bethe ansatz solutions and critical properties. As far as we 
know, the only exception is the UqOsp{2\2) |^ model which has been recently solved in 
the context of = 2 supersymmetric lattice models 0. 

The purpose of this letter is to exactly solve the simplest system in the Osp{n\2'm) 
family, namely the case of one boson {n = 1) and two fermions (m = 1). As we shall see, 
a new feature appears in the Bethe ansatz equations, making possible an explicit distinc- 
tion between the even and odd sectors of the model. This new property has important 
consequences for the finite-size behaviour which governs the class of universality of this 
theory. 

We begin by introducing the isotropic Osp{l\2) magnet. Its Hamiltonian defined on a 
lattice of L sites (assuming periodic boundary conditions) can be written in the following 
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compact form 



^=-E[^M+l + TT^M+l] 
i=l 



where {Pi,i+i)f 



lab — ^a,dSb,c is the graded permutation operator and {Ei^i+i)lf = 

aabO^cd is the Osp{l\2) Temperely-Lieb invariant operator ]10[. The variables p{a) are the 
Grasmann parities, namely p{l) = (boson) and p{2) = p{3) = 1 (fermions). The symbol 
a** denotes the supertranspose operation on a 3X3 matrix a given by 
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The Bethe ansatz solution of (1) is as follows. The first step is to write the cor- 
responding vertex operator of the lattice model leading to the Hamiltonian (1). After 
an appropriate canonical transformation, it is possible to show that this weight can be 
expressed in terms of the generators S*^, ^ S~ of a spin-1 algebra. Its explicit form is 



R{\) 



^[S'S+ ~2Xh{X)S+S'] [X - 3h{X)]I + 3h{X){S'f -^[S'S' + 2Xh{X)S-S' 

/(A)(S+)2 -j^[S^S^ + 2XhiX)S^S+] XI - fiX)S^ - g{X)iS'y 

(3) 

where / is the 3X3 identity matrix and functions /(A), g{X) and h{X) are given by 



V 



/(A) - MA) 



2(xhj2y ^(^)-2A + ^MA) 



(4) 



2(A-3/2) 

It is important to notice that the reference state preserving the triangularity property of 
this i?-matrix is still the usual ferromagnetic pseudo-vaccum. Moreover, we recall that its 



structure resembles much the one that appears in the Izergin-Korepin (IK) model |jTT], . 
Indeed we are able to show that the Hamiltonian (1) can be obtained as a certain branch 
limit of the IK system through a delicate but rigorous canonical transformation. The 
calculations leading to this result are rather cumbersome and will be presented elsewhere 



[O]. Here we can use this information ,by adapting the algebraic and(or) the analytic 



Bethe ansatz method used in the IK model |jT2|, in order to solve Hamiltonian (1). Due to 
a f/(l)-conserved charge the Osp{l\2) Hamiltonian can be block diagonalized into disjoint 
sectors labelled by the total magnetization index, namely r = Ylf=i S^. In a given sector 
r, for periodic boundary condition, the eigenstates are parametrized by a set of complex 
numbers Xj satisfying the following Bethe ansatz equation 

and the eigenenergies are given by 

The important novelty present in the Bethe ansatz equation is its dependence on a 
phase factor (— 1)*^ even for a system with periodic boundary conditions. Hence, even and 
odd sectors of the theory are distinguished by the signs of their "bare" phase-shift (right 
hand side of (4) ). This may be interpreted as an explicit separation of the bosonic and 
fermionic degrees of freedom present in the system. In fact, we have verified this "sector 
separation" through an exact diagonalization of the Hamiltonian (1) for small values of L. 
We also remark that although a phase factor does not affect the basic properties in the 
thermodynamic limit {L oo), they do change the finite-size behaviour. Therefore, this 
phase factor plays an important role in order to characterize the correct critical behaviour 
of the model. 

Before turning to the study of the finite size properties of (4) and (5), we first compute 
some useful quantities in the thermodynamic limit. The ground state in each sector r 
is parametrized by a set of real numbers A^. Taking the logarithm of equation (4) and 
correctly collecting the phase factors, we rewrite the Bethe ansatz equations as 

L-r 

L^i/2{Xj) = 27rQj + lM>^j - >^k) - V'i/2(A, - A^)] (7) 

where ipa{x) = 2arctan(x/a) and Qj are integer or semi- integer series of number 

Qj = - '^~^~^^ J = l,2,---,L-r (8) 
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For large L, the roots \j are densely packed into its density distribution, crii^j) — 
l/L(Aj+i — Strictly in the L — oo limit this system can be solved by elementary 
Fourier techniques and we find that 

^ 2 cosh(27rA/3) 
^ ' v^cosh(47rA/3) + l/2 ^ ^ 

The ground state energy (per site) can be calculated by using equations (6) and (9) 
after replacing the sum by an integral. The final result is 

Coo = lim Eo{L)/L = — ^ + 1 = -1.4184 (10) 

L^oo 9 

In addition, we have verified that the low-lying excitations over the ground state are 
gapless. As usual, the excitations are built up by inserting vacancies on the density 
distribution of Xj. The simplest excitation is the triplet state made of two real "holes" Ai 
and A2. Their contribution to the total momenta is calculated to be 

2 

p{Xi, X2) = TT — 2 arctan[sinh(27rAi/3)/ cosh(7r/6)] (11) 

4 = 1 

which leads to the following low-momentum dispersion relation 

e(p) ~ VsP (12) 

where = 27r/3 is the sound velocity. 

At this point we have obtained the most important quantities necessary for the analysis 
of the critical properties. The class of universality can be determined by exploiting a set of 
important relations between the eigenspectrum of finite-lattice system [jl^. For instance. 



the central charge c is related to the ground state energies Eq{L) by [ P^ 

Eo{L)/L-e^ = -^ (13) 

In Table 1 we present our results for the central charge c. This result is obtained by 
numerically solving the Bethe ansatz equation in the sector r = 0. The extrapoled results 
predict a conformal anomaly c = 1. The sectors r 7^ are responsible for the low- lying 



conformal dimensions. In this case , however, a numerical study up to lattice size L ~ 64 
presents a poor convergence rate. This is due to the appearance of logarithm corrections 
of order 0(l/ln(L)) in the scaling function. Fortunately, since the roots Xj are real, we 
can use an analytical method developed by De Vega and Woynarovich [jl6l in order to 
estimate these critical dimensions. The calculation is quite technical |jl3[ and here we only 
present the final results. We stress, however, that the phase factor of (5) is of fundamental 
importance in this computation. The critical dimension Xr associated with the lowest 
operator on sector r has the following form 

Xr = - (14) 

Such dimension corresponds to a spin- wave excitation of index r in a Gaussian structure 
of operators. In fact, more generally we expect to find the following dimensions 

Xr,m = —+m^ (15) 

where m stands for the "vortex " excitations index. On the other hand, this result can 
be justified by using the following arguments of symmetry. For example, we have verifed 
that ( for finite lattices) the first excited state with zero momentum in the r = sector 
is degenerated to the ground state in sector r = 2. This means the identity Xo,i = X2fl, 
which is indeed verified in formula (15). Evidently, the same argument can be repeated 
for many other values of r and m, confirming the conjecture (15). We recall that this is 
the manifestation of the Osp{l\2) symmetry encoded in the important phase factor (— 1)^. 

In summary, we have given a Bethe ansatz solution to the eigenvalues of the Osp{l\2) 
magnet. The novelty is the presence of a phase factor depending on the sector r, which 
becomes crucial in order to determine its critical properties. We believe that our results 
extends the form of Bethe ansatz solution appearing in integrable lattice models. Hope- 
fully, our results will be of relevance in the discussion of the Bethe ansatz completeness 
and in the thermodynamic properties. 



5 



Acknowledgements 

It is a pleasure to thank M. Malvezzi for his help with numerical checks and F.C. Alcaraz 
for innumerous discussions. This work is supported by CNPq (Brazilian agency). 

References 

[1] p.p. Kulish and E.K. Sklyanin, J.Sov.Math.19 (1982) 1596 

[2] V.V. Bazhanov and A.G. Shadrikov, Theor.Math.Phys.73 (1987) 1302 

[3] J.H.H. Perk and C.L. Schultz, Non-Linear Integrable Systems- Classical and Quantum 
Theory, ed M. Jimbo and T. Miwa (World Scientific, Singapore) pl35 

[4] H.J. de Vega and E. Lopes, Phys.Rev.Lett.67 (1991) 489 

[5] P. Schlottmann, Phys.Rev.B 36 (1987) 5177 

P.A. Bares and G. Blatter, Phys. Rev. Lett. 64 (1990) 2567 

[6] F.H.L. Esslcr and V.E. Korepin, Phys.Rev.B 46 (1992) 9147 

F.H.Essler, V.E. Korepin and K. Schoutens, Phys.Rev.Lett.68 (1992) 2960, 
Phys. Rev. Lett. 70 (1993) 73 , and references therein 

[7] A. Foerster and M. Karowski, Nucl.Phys.B 396 (1993) 611 

[8] T. Deguchi, A. Fujii and K. Ito, Phys.Lett.B 238 (1990) 242 

[9] Z. Maassarani, preprint-USC-94/007, april 1994 

[10] M.J. Martins and P.B. Ramos, preprint-UFSCARF-TH-94-12, to be published in 
J.Phys.A 

[11] A.G. Izergin and V.E. Korepin, Commun. Math. Phys. 79 (1981) 303 

6 



[12] V.I. Virchirko and N.Yu. Reshetikhin, Teor.Mat.Fiz. 56 (1983) 260 
V.A. Tarasov, Theor.Mat.Fiz. 76 (1988) 793 

[13] M.J. Martins, in preparation 

[14] J.L. Cardy, Phase Transitions and Critical Phenomena vol. 11, ed C. Domb and J.L. 
Lebowitz (New York: Academic) p55 

[15] H. B16te, J.L. Cardy and M.P. Nightingale, Phys.Rev.Lett.56 (1986) 742 
I. Affleck, Phys.Rev.Lett. 56 (1996) 746 

[16] H.J. de Vega and F. Woynarovich, Nucl.Phys.B 251 (1985) 439 

Table Captions 

Table 1. The estimative of the conformal anomaly of equation (13) 



Table 1 



L 


-6L^[Eo{L)/L-e^]/nVs 


16 


1.0029 729 


24 


1.0022 324 


32 


1.0018 406 


40 


1.0015 955 


48 


1.0014 261 


56 


1.0013 010 


64 


1.0012 043 


Extrapolated 


1.0001 (1) 
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